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Abstract 

We develop a general framework, based on Boltzmann transport 
theory, to analyze the distribution of wealth in societies. Within this 
framework we derive the distribution function of wealth by using a 
two-party trading model for the poor people while for the rich peo- 
ple a new model is proposed where interaction with wealthy enti- 
ties (huge reservoir) is relevant. At equilibrium, the interaction with 
wealthy entities gives a power-law (Pareto-like) behavior in the wealth 
distribution while the two-party interaction gives a Boltzmann-Gibbs 
distribution. 
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1 INTRODUCTION 

Inequality in the distribution of wealth in the population of a nation has 
provoked a lot of studies. It is important for both economists and physicists 
to understand the root cause on this inequality: whether stochasticity or a 
loaded dice is the main culprit for such a lop-sided distribution. While it 
has been empirically observed by Pareto |]Q that the higher wealth group 
distribution has a power-law tail with exponent varying between 2 and 3, 
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the lower wealth group distribution is exponential or Boltzmann-Gibb's like 
[3 IS]- The Boltzmann-Gibb's law has been shown to be obtainable when 
trading between two people, in the absence of any savings, is totally random 
jH 03 |B] . The constant finite savings case has been studied earlier numeri- 
cally by Chakraborti and Chakrabarti 0j and later analytically by us [Hj- As 
regards the fat tail in the wealth distribution, several researchers have ob- 
tained Pareto-like behavior using approaches such as random savings |S] , 
inelastic scattering 0|, generalized Lotka Volterra dynamics JT0], analogy 
with directed polymers in random media [TT], and three parameter based 
trade-investment model [T2] . 

In this paper, we try to identify the processes that lead to the wealth 
distribution in societies. Our model involves two types of trading processes 
- tiny and gross. The tiny process involves trading between two individuals 
while the gross one involves trading between an individual and the gross- 
system. The philosophy is that small wealth distribution is governed by 
two-party trading while the large wealth distribution involves big players 
interacting with the gross-system. The poor are mainly involved in trading 
with other poor individuals. Whereas the big players mainly interact with 
large entities/organizations such as government (s), markets of nations, etc. 
These large entities/organizations are treated as making up the gross-system 
in our model. The gross-system is thus a huge reservoir of wealth. Hence, 
our model invokes the tiny channel at small wealths while at large wealths 
the gross channel gets turned on. 

2 GENERAL FRAMEWORK 

We will now develop a formalism similar to Boltzmann transport theory 
so as to obtain the distribution function f(y,y,t) for wealth y, net income 
y (or total income after consumption) as a function of time t. Similar to 
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Boltzmann's postulate we also postulate a dynamic law of the form 

V. J ext source < J interaction 

The first term on the right hand side (RHS) describes the evolution due to 
external income sources only, while the second one represents contribution 
from entirely internal interactions. 

2.1 Model for tiny-trading 

Individuals, possessing wealth smaller than a cutoff wealth y c , engage in 
two-party trading where two individuals 1 and 2 put forth a fraction of their 
wealth (1 — \t)yi and (1 — Xt)y2 respectively [with < \ t < 1]. Then the 
total money (1 — \t)(yi + 2/2) is randomly distributed between the two. The 
total money between the two is conserved in the two-party trading process. 
We assume that probability of trading by individuals having certain money 
is proportional to the number of individuals with that money. 

Since in a two-party trading there is no external income source, then 
{ % \ =0. The second term on the RHS in Eq. (JTJ) can be obtained 

L at J ext source 

as follows in terms of a balance equation. 



?L 

dt 




= gains — losses. (2) 

interaction 



In Eq. (J2J) the two terms on the RHS can be expressed in terms of transition 
rates r(yi, y 2 ; y[, y' 2 ) for a pair of persons to go from moneys y±, yi to moneys 
y[, y' 2 respectively. Then, we have on assuming that the distribution function 
is only a function of wealth and time, 

0/(3/i,*) 



dt 



r(y[, y' 2 ] Vi, 2/2)/(2/i, t)f(y' 2 , t)dy 2 dy[dy 2 
r(yu 2/2; y[, y' 2 )f(yi, t)f(y 2 , t)dy 2 dy[dy' 2 . (3) 



In the above equation, conservation law requires that 2/1 + 2/2 = 2/1 + 2/2- Hence 
in the first integral we treat y 2 as redundant and integrate out with respect 
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to it to yield a normalization constant. Similarly in the second integral y' 2 is 
integrated out. Now for the transition rate in the first integral of the above 
equation we have 

r(s! '* s " rio: (i-W 1+ *' (4) 

if Ky[ < Vi < y[ + (1 — ^t)y ! 2 an d zero otherwise. On taking into account the 
restriction that no one can have negative money and setting y[+y' 2 = L, the 
first integral in Eq. (J3J) at equilibrium is proportional to 

rco rb(y 1 ,L,\ t ) 

/ dL dy' 1 J r (y' 1 , L, X t ), (5) 

Jyi Ja(yi,L,X t ) 

where 

a(yi,L,X t ) = max[0,{yi - (1 - \ t )L}/X t ], 



6(2/1, L, A*) = mm[L,y 1 /X t ], 

and 

m)f(L-y{) 



X f )L 



As regards the second integral in Eq. (jSJ), at equilibrium we assume that the 
transition from y\ to all other levels is proportional to f(yi)- Also, since at 
equilibrium d ^Q^^ = 0, we obtain the distribution function to be 

roc rb(y,L,X t ) 

f(y)= / dL / dxJ 7 (x,L,X t ). (6) 

The above result was obtained earlier by us using an alternate route [3]. 

On introducing an upper cutoff y c for the two-party trading, the contri- 
bution to the distribution function f(y) from the tiny channel becomes 

roo My,L,X t ) 

7 / dL dxJ 7 (x : L,X t )H(x,L,y c ). (7) 

Jy Ja(y,L,\ t ) 
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In the above equation 

H(x, L, y c ) = [1 - 9(x - y c )\ [1 — 6(L — x — y c )}, 

with 6(x) being the unit step function and 7 = 1/ f$ c dxf(x) is a normal- 
ization constant introduced to account for the less than unity value of the 
probability of picking a person below y c . 

2.2 Model for gross-trading 

Next, we will analyze the contribution to the distribution function f(y) from 
gross-trading. An individual possessing wealth larger than a cutoff wealth 
(y c ) trades with a fraction (1 — X g ) of his wealth y with the gross-system. The 
latter puts forth an equal amount of money (1 — X g )y. The trading involves 
the total sum 2(1 — \ g )y being randomly distributed between the individual 
and the reservoir. Thus on an average the gross-system's wealth is conserved. 
When only the gross-trading channel is operative, we have 

dJ ^- = Jr(y[;y 1 )f(y[,t)dy' 1 

-J r(yi;y[)f( yi ,t)dy[, (8) 

where r(y[; y{) is the transition rate from y[ to y\ through interaction with 
the reservoir. Total money involved in trading (between individual and the 
gross-system) is 2^(1 — \ g ). After interaction, the resulting money y\ of the 
individual satisfies the following constraints X g y[ < y\ < (2 — X g )y[. The 
first transition rate in the above equation is given by 

if \ g y[ < yi < (2 — \ g )y[ and zero otherwise. As before, at equilibrium 
the second integral in Eq. (JHJ) is proportional to f{yi) when only the gross 
channel is operative. Then the distribution function f(y) is given by 

= rv/* B dxf(x) 
J{y) Jy/(2-x g ) 2x(l - X g y { ' 
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Now it is interesting to note that the solution of the above equation is given 
by f{y) = c/y n . To obtain n one then solves the equation 

(2-X 9 r-X n g =2n(l-X g ), (11) 

and obtains n = 1,2. Only n = 2 is a realistic solution because it gives a 
finite cumulative probability. It is of interest to note that the solution is 
independent of X g . Also, clearly the distribution function makes sense only 
for y > 0. On taking into account an upper cutoff y c , the contribution to the 
distribution function f(y) from the gross channel is 

rv/ x 9 dxf(x) n . . 

9(x-y c ). (12) 



ly/(2-X g ) 2x(l - X g ) 

2.3 Hybrid model 

Here an individual possessing wealth larger than a cutoff wealth y c does 
trading with the gross-system, while individuals possessing wealth smaller 
than y c engage in two-party tiny-trading. Hence from Eqs. (JJJ) and ()12J) . the 
distribution function is obtained to be 

foo rb(y,L,\ t ) 

f(y)= 7/ dL dxJ r (x, L, X t )H(x, L, y c ) 

rvA» dxf(x) 

Jy/(2~\ g ) 2x(l - Xg) 

Now, it must be pointed out that when the savings X t = 0, X g ^ 0, and 
y — > 0, Eq. (|T3j) yields (up to a proportionality constant) the following same 
result as the purely tiny-trading case without an upper cutoff [5j: 

f'(y) cc -f(y)f(0). (14) 

In obtaining the above equation we again assumed that the function f(y) 
and its first and second derivatives are well behaved. Then the solution for 
small y is given by 

f(y) <x f(0)exp[-yf(0)]. (15) 
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3 RESULTS AND DISCUSSION 

The distribution function f(y) can be obtained by solving the nonlinear inte- 
gral Eq. (|13p. To this end, we simplify Eq. (JTSJl for computational purposes 
as follows: 

r 2y c My,L,M) 



rzy c rO(y,^,A t ) 

f(y)= lG{y,h,yc) \ dL dxJ 7 (x,L,X t )H(x,L,y c ) 

Jy Ja(y,L,\t) 
r-. n, m f y/Xs dxf(x) „, . 

+ 1 - % - y os / - ; x - y c ) 

Jy/(2-X g ) 2x(l - X g ) 



+0(y-y as )f(yas)^i, (16) 
y 

where C/(y, X t ,y c ) = 1 — 0[y — (2 - \ t )y c } and y > y as gives the asymptotic 
behavior /(y) oc 1/y 2 . In our calculations, we have taken y as to be at least 
20y c and obtained f(y) for all y less than 2000 times the average wealth 
per person y av . We solved Eq. (fTT)|) iteratively by choosing a trial function, 
substituting it on the RHS and obtaining a new trial function and successively 
substituting the new trial functions over and over again on the RHS until 
convergence is achieved. The criterion for convergence was that the difference 
between the new trial function /„ and the previous trial function f p satisfies 
the accuracy test \f n {Vi) ~ fp{Vi)\/T,i fp{Vi) < 0.002 [I3J. 

In Fig. 1, using a log-log plot we depict the distribution function f(y) for 
the constant savings case X t = \ g = 0.5 with the average money per person 
y av being set to unity and with the values of the wealth cutoff y c = 3, 5, 10. 
As expected, for larger values of y c , the Pareto-like 1/y 2 behavior sets in later. 
The transition to purely gross-trading occurs at (2 — Xt)y c , while below X g y c 
it is purely two-party tiny-trading. Thus the transition from purely tiny- 
trading to purely gross-trading occurs in Fig. 1 over a region of width y c . 
However, all the tails merge irrespective of the cutoff. At smaller values of y 
the behavior of f(y), depicted in the inset, is similar to the purely two-party 
trading model studied earlier (see Ref. 0). The curves in the inset appear 
to be close because here the trading is two-party and is governed by the same 
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Figure 1: Plot of the wealth distribution function for savings At = X g = 0.5 
and various wealth cutoff values y c — 3, 5, 10. The average money per person 
y av is set to unity. The dotted lines are guides to the eye. 
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Figure 2: Wealth distribution f(y) at average wealth y av = 
y c = 5, and various values of savings A t = \ g = 0.1, 0.5, 0.8. 



1, wealth cutoff 
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Figure 3: Money distribution function at zero savings for tiny-trading and 
various savings values X g = 0.2, 0.5, 0.9 for the gross-trading. The average 
money y av = 1 and the wealth cutoff y c — 5. 

savings. Next, in Fig. 2 we plot f(y) with the cutoff y c = 5, y av = 1, and 
for values of savings fraction A t = X g = A = 0.1, 0.5, 0.8. Here the power-law 
behavior {1/y 2 ) takes over for y > (2 — A)?/ c and hence at lower savings it 
sets in later. In the power-law region the curves merge together. As shown 
in the inset of Fig. 2, at smaller values of y the f(y)s become zero with the 
higher peaked curves (corresponding to larger As) approaching zero faster 
similar to the case of the purely two-party trading model in our earlier work 
Here the transition from purely tiny- to purely gross-trading at higher 
A is sharper because the transition occurs over a region of width 2(1 — X)y c . 
Lastly, in Fig. 3, we show the distribution function f(y) for the zero savings 
case in the tiny-channel (Xt = 0) and for various savings X g = 0.2, 0.5, 0.9 in 
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the gross-channel with y av = 1 and y c = 5. The distribution, as expected, 
decays exponentially (or Boltzmann-Gibbs-like) for small values of y and has 
power-law (1/y 2 ) behavior at large values. The curves merge in the Pareto- 
like region and, in fact, f(y) ~ 0.1/y 2 in all the three figures at large values 
of y. In Fig. 3 too, for reasons mentioned earlier, the transition is sharper at 
larger values of A g . Fig. 3 takes into account the fact that, in societies, the 
rich tend to have higher savings fraction (A) compared to the poor. Actually, 
if the savings fraction were to increase gradually with wealth, one can expect 
a more gradual change in the transition region of the distribution rather than 
the sharp local maxima (around y « 6.5 ) shown by the X g = 0.9 curve. 

In all the figures anomalous looking kinks/shoulders appear at the cross 
over between the Boltzmann-Gibbs-like and the Pareto-like regimes. This 
is due to the sharp cut-off at y c that we introduced using a step function. 
However, a kink seems to be generic in these kind of distributions in real pop- 
ulations (as borne out by the empirical data in Fig. 9 of Ref. indicating 
that two different dynamics may be operative in the two regimes. Different 
societies have the onset of Pareto-like behavior at different wealths which is 
indicative that the cut-off has to be obtained empirically based on various 
factors like the social structure, welfare policies, type of markets, form of 
government, etc. 

In Japan the wealth/income distribution vanishes at zero wealth/income 
and then rises to a maximum (see Ref. [21). In US the distribution seems to 
be a maximum at zero wealth/income (see Ref. j2j). Both these aspects can 
be covered in our model as the poor in general are known to save very little. 
If their savings are zero, one gets the Boltzmann-Gibbs behavior at the poor 
end. On the other hand, if the savings are small one gets a maximum close 
to zero and the distribution vanishes at zero wealth. 

It would be interesting to deduce the savings pattern from the wealth 
distribution. While it has been observed that the rich tend to save more 
than the poor, how gradually the savings change as wealth increases can 
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perhaps be inferred from the change in slope. However, as explained below, 
the middle region (involving the middle-class) has been modeled quite crudely 
by us and needs to be refined before a serious connection with savings pattern 
can be attempted. 

We will now further discuss the motivation for using two different mech- 
anisms to model the observed wealth distribution. The model is an approx- 
imation where the direct wealth exchange occurs between people who are 
in economic proximity. At the poorer end of the spectrum, the poor, who 
have limited economic means and avenues, come in contact with a few poor 
and their economic activity is modeled in terms of two-party trading. At 
the other end of the wealth spectrum, the rich have access to various eco- 
nomic avenues (e.g., markets, know-how, work force, capital, credit facilities, 
contacts, wealthy society, etc.) due to which they can trade with huge or- 
ganizations and are thus modeled to interact with a reservoir. As regards 
the middle-class that is between the rich and the poor, they trade amongst 
themselves as well as with the poor and the reservoir. As a first step towards 
realizing this scenario, we included in our model only the two extreme cases 
of interaction. What we have not taken into account is the interaction of 
the middle class with the reservoir. To rectify this, in future we hope to 
introduce a cutoff y g for the interaction with the reservoir such that y g lies 
below the two-party trading cutoff y t . Thus, we believe that our model is a 
reasonable one at the poor and rich ends and is a crude approximation for 
the middle class. While it is true that the poor also come in market contact 
with wealthy organizations like the coke company, the contact is an indirect 
one through intermediaries. For example, the poor person deals with a richer 
shop-keeper selling coke who in turn deals with a richer local distributor who 
in turn deals with the big coke company. Lastly, we would like to add that 
the assumption of random distribution in two-party trading is a model stud- 
ied by others as well (see Refs. [HE]). We feel that in any trading there is 
random fluctuation of the price around its true value. The total money put 



13 



forth for trading corresponds to the amount of random fluctuation. However 
the poorer of the two puts forth less and makes the trading biased in his/her 
favor. This can be justified from the fact that the poor people are constantly 
looking for bargains to make ends meet. 

Compared to other types of analyses involving two-party trading to ex- 
plain Pareto law (see Refs. [HE]), our g ross -trading mechanism can make 
contact with the standard approach in macroeconomics as will be shown 
below. Over the past, economists have developed two models, namely, the 
dynastic model and the life-cycle model, to explain wealth distribution. In 
the dynastic model, where bequests are vehicles of transmission of wealth 
inequality, people save to improve the consumption of their descendants. On 
the other hand, in the life-cycle model, where wealth of an individual is a 
function of the age, people save to provide for their own consumption after re- 
tirement. Both these models and their hybrid versions have had only limited 
success |14j . However, one of the ingredients that goes into these models, 
i.e., uninsurable shocks or stochasticity in income, has been exploited by 
econophysicists with remarkable success in reproducing power-law tails. 

In macroeconomics, the objective is to maximize a cumulative utility func- 
tion subject to a wealth constraint ^H]. Mathematically this is formulated 
as 

max E^PHct+i), (17) 
ct+i,yt+i '—' 
i 

subject to the constraint 

y t+i = (1 + r)y t+ j_i + e t+i - c t+h (18) 

where c t , yt, and e t are consumption, wealth, and labor earnings respectively 
at time t, r is the interest rate on wealth y, < (3 < 1 is the time-discount 
factor, u(ct) is the concave utility function, E t is the expectation value based 
on the available information at time t. Using the method of Lagrange mul- 
tipliers, the conditions of optimality yield 

E t [u\c t ) - (1 + OMct+i)] = 0, (19) 
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where u'(ct) is the derivative of u(ct) with respect to q. From the above 
equation we see that consumption at different times are related. In our work 
[see Eq. Q], we introduced the stochasticity 

Vt+i ~Vt = e (! - h)Vti (20) 

where e is a random number with — 1 < e < 1, which implies that 

ry t + e t+1 - c t+1 = e(l - X g )y t . (21) 

The above equation can be made consistent with the optimal consumption 
relation given by Eq. (|T9"|) . Thus our results can be approached through the 
standard machinery in macroeconomics. 

The stochasticity in wealth given by Eq. ()20|) implies that the spread 
in wealth distribution at time t + 1 is proportional to y t and thus wealth's 
y t +i > yt yield a wider spread for y t+2 than do wealth's y t +i < yt- Thus the 
distribution becomes more skewed to the right. 

In conclusion, we introduced a new ingredient - interaction of the rich 
with huge entities - and obtained a Pareto-like power-law. On the other 
hand, the Boltzmann-Gibbs-like wealth distribution of the poorer part of the 
society is explained through a two-party trading mechanism. All in all, we 
demonstrate that stochasticity can account for the observed skewness in the 
wealth distribution. 
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